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The relative velocity correlation function of pairs of intermediate mass fragments has 
been studied for d+Au collisions at 4.4 GeV. Experimental correlation functions are 
compared to that obtained by multibody Coulomb trajectory calculations under the 
assumption of various decay times of the fragmenting system. The combined approach 
with the empirically modified intranuclear cascade code followed by the statistical multi- 
fragmentation model was used to generate the starting conditions for these calculations. 
The fragment emission time is found to be less than 40 fm/c. 
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1. Introduction 

The main decay mode of very excited nuclei (E* > 4 MeV/nucleon) is copious emis- 
sion of intermediate mass fragments (IMF), which are heavier than a-particles but 
hghter than fission fragments. An effective way to produce hot nuclei is reactions 
induced by heavy ions with energies up to hundreds of MeV per nucleon. But in this 
case the heating of the nuclei may be accompanied by compression, rotation, and 
shape distortion, which can essentially influence the decay properties of hot nuclei. 
The picture becomes clearer when light relativistic projectiles (protons, antipro- 
tons, pious) are used. In this case, fragments are emitted by only one source (the 
slowly moving target spectator). Its excitation energy is almost entirely thermal. 
Light relativistic projectiles provide therefore a unique possibility for investigating 
thermal multifragmcntation. 

It has been found experimentally that the process is characterized by two vol- 
umes ^i^. The first one corresponds to the chemical freeze-out state when prefrag- 
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ments are formed, Vt = (2.6 ± 0.2)Vb. The second one is reached by the nucleus 
after its descent from the top of the fragmentation barrier to the multi-scission 
point. It is cahed the kinetic freeze-out volume Vf = (5.0 ± 0.5)Vb. 

The decay pro perti es of hot nuclei are weU described by statistical models of 
multifragmentatiorPE! and this can be considered as an indication that the system 
is in thermal equilibrium or at least close to that. 

The time scale of fragment emission is a key parameter for understanding the 
decay mechanism of highly excited nuclei. Is it sequential and independent evapo- 
ration of IMFs or is it a multibody decay mode with almost simultaneous emission 
of fragments governed by the total accessible phase space? As was suggested by 
D.H.E. Gross in ref.^^ "simultaneous" means that fragments are liberated during 
a time interval which is smaller than the Coulomb acceleration time Tc, when the 
kinetic energy of fragments amounts to 90% of the asymptotical value. According 
to^^, Tc is estimated to be (400 - 500) fm/c. Fragments emitted within this time 
interval are considered being not independent as they interact via the Coulomb 
force while being accelerated in the electric field of the source. As a result, the 
yield of events with small relative velocities of the fragments (or small relative an- 
gle between them) is suppressed. The magnitude of this effect drastically depends 
on the emission time since the longer the time separation of the fragments, the 
larger their space separation and the weaker the mutual Coulomb repulsion. Thus, 
measurement of the IMF emission time Tem (the mean time separation between two 
fragment emissions in a given event) is a direct way to answer the question as to the 
nature of the multifragmentation phenomenon. In some papers, the mean lifetime 
of the fragmenting system is used to characterize t he p rocess of disintegration. 
There is a simple relation between these two quantities''^ 



Both values are close to each other when the mean IMF multiplicity M is in the 
range 2-3, as in the case of the light relativistic projectiles. 

The time scale for IMF emission is estimated by comparing the measured cor- 
relation function with the multibody Coulomb trajectory calculations with Tem as 
a parameter. There are two procedures to measure the emission time: analysis of 
the IMF-IMF correlation function of the relative angle or the relative velocity. 

The first measurements of the time scale for the thermal multifragmentation 
were performed for ^He+Au collisions at 14.6 GeV by analyzing the IMF-IMF 
relative angle correlatiorP^^. It was found that Tem is less than 75 fm/c. The same 
procedure was used by the FAS A group for the p+Au interaction at 8.1 
when the emission time Tem < 70 fm/c was found. A similar value was obtained 
by the ISiS collaboration for collisions of 4.8 GeV "^He with a gold targe P. In this 
paper IMF-IMF relative velocity correlations were studied. A general overview of 
the experimental activity in this field can be found in review paper ^1^. 

In the present paper the temporal characteristics of multifragmentation are 
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investigated for the first time for interaction of 4.4 GeV deuterons with the Au 
target. For that purpose the IMF-IMF correlation function of the relative velocity 
was studied. 

2. Experimental 

The experiment has been performed with the Att setup installed at 

the external beam of the Dubna superconducting accelerator NUCLOTRON. The 
FASA device consists of two main parts: 

• the array of thirty AE-E telescopes, which serve as triggers for the read-out of 
the whole FASA detector system. These telescopes allow measuring the fragment 
charge and energy. The spatial distribution of fragments is also obtained (with 
steps Ae = 10°), 

• the fragment multiplicity detector (FMD) including 58 thin CsI(Tl) counters 
(with scintillator thickness around 35mg/cm^), which cover 81% of 47r. The FMD 
gives the number of IMF's in the event and their spatial distribution. 

The fragment telescopes consist of a compact ionization chamber as the AE 
counter and a Si (Au) semiconductor detector as the E spectrometer. Effective 
thickness of the E detector was around 700/i, which is enough to measure the energy 
spectra of all intermediate mass fragments. The ionization chambers have a shape 
of a cylinder (50 mm in diameter, 40 mm in height) and are made from polished 
brass. The entrance and exit windows are made from organic films (^100/ig/cm^) 
covered by a thin gold layer prepared by thermal evaporation. A gold wire 0.5 mm 
in diameter is used as the anode. The cathode (brass cylinder and mechanically 
supported thin entrance and exit windows) is surely grounded. Carbon fluoride 
CF4 at the pressure of 50 torr is used as a working gas. 

A self-supporting Au target (^1.5 mg/cm^) is located at the center of the FASA 
vacuum chamber supported by thin tungsten wires. The energy calibration of the 
counters was done periodically using a precise pulse generator and a ^"^^Am a 
source. The beam intensity was around 10^*^ particles per spill. The beam spot 
was continuously controlled by two multiwire proportional chambers placed at the 
entrance and the exit of the FASA device. The beam intensity was measured by 
the special ionization chamber located 150 cm behind the target. The spill length 
was 1.5 s, the frequency of the beam bursts was 0.1 Hz. 

Fig. 1 gives an example of the AE-E plot measured by one of the telescopes for 
d-|-Au collisions at the beam energy 4.4 GeV. The energy is given for the fragment 
at the front of the telescope. It was obtained from the measured value by adding the 
energy absorbed on the way from the target to the E-detector. The loci for different 
intermediate mass fragments are well resolved in the range from lithium to silicon. 
The energy cutoff caused by adsorption of the low energy fragments before the 
Si(Au) detector increases with Z of fragments. Thus, Ne fragments with the energy 
lower than 25 MeV do not reach the E-counter. 
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Fig. 1. AE-E plot measured by one of the thirty telescopes of FASA device, total number of events 
detected is around 1.5- 10®. 



3. Results of Measurements 

We used a refined version of the intranuclear cascade model (INC to describe the 
fast stage of the reaction and to get the N, Z and the excitation energy distributions 
of the target spectators. The mean exc itation energy of the target spectator is 
(150-200) MeV according to our analysiJI^H^. The excitation energy distribution 
is rather wide, and 10-15% of the total inelastic cross-section correspond to the 
excitation energies of residuals higher than the fragmentation thr esho ld (300-400 
MeV). We use the Moscow- Copenhagen Statistical Model (SMm|SI11 to describe 
the disintegration of hot residuals. The primary fragments are hot, and their de- 
excitation is also considered by the SMM to get the final distributions of cold 
IMF's. The last stage of the reaction is Coulomb expansion of the final system, 
when fragments are getting their kinetic energy. This process is described by the 
multibody trajectory calculations for fragments in the freeze-out volume Vf — 3Vq. 
This part of the combined model gives the IMF kinetic energy spectra and velocity 
distributions for fragments. 

It was shown in our paperJ^IIll that traditional (INC + SMM) model calcu- 
lations failed to describe the data for the IMF multiplicities. It was considered as 
anindication that the cascade calculations overestimate the excitation energies of 
the residues. In order to overcome this difficulty the calculated excitation energies 
are reduced empirically by the factor a on the event-by-event basis 



a = — — ^ (2) 

< MiNC+SMM > 

where < Mexp > - measured mean multiplicities for events with at least one IMF. 
The model calculated mean IMF multiplicity is given as a denominator. The mass 
loss during the intra-nuclear cascade is reduced by the factoif^^ (1 - a). As a result, 
one gets a INC-|-Exp-)-SMM model which we used in the present paper. 
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Fig. 2. The yield as a function of the reduced relative velocity of coincident IMF's (points). 
Histogram is for the mixed events. 



Figure 2 shows the distribution of the reduced relative velocity of the fragments. 

Vred = Vrel/V Zi + Z2 (3) 

Difference of velocities of two fragments from the same event is normalized by the 
square ro ot of the sum the fragment charges. This normalization was suggested in 
| 16 | 17 | 18 | compensate charge dependence of the velocity difference | Vi-Vk |. It is 
measured for 18129 events with detection of at least two coincident fragments (the 
total number of triggers was 4-10^). Data for all the fragments with Z between 3 
and 20 are included. 

The yield of a pairs with small Vred is remarkably suppressed. It is partly related 
to the phase space factor (or phase space density), as is demonstrated in fig. 2 by 
the histogram obtained for the "mixed events" (i.e. for a pair of fragments from 
different collisions). But the Coulomb repulsion between the coincident fragments 
makes suppression at Vred ~ much more significant. 

Traditionally, Coulomb interaction of genetically related fragments is described 
in the terms of the correlation function 



1 + R{Vred) — C Nco7-r{Vred) / ]^uncorr{Vred) : 



(4) 



where Ncorr{vred) IS the measured coincidence yield. The denominator is composed 
of the mixed events, C is the normalization constant. 



4. Discussion 

The model correlation functions are obtained after multibody trajectory calcula- 
tions for the fragments in the Coulomb field of the system. It is made on an event by 
event basis using the INC+Exp+SMM model for describing specification of events 
(Z, A and excitation energy of IMF source). 
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Fig. 3. Reduced relative velocity correlation function of the coincident fragments pairs. The points 
are experimental data, lines are obtained with INC+Exp+SMM combined model followed by the 
multi-body Coulomb trajectory calculations with decay time of fragmenting system: a - prompt, 
b - 160 fm/c. 

The starting time for each fragment to move along a Coulomb trajectory was 
randomly chosen according to the decay probability of the system, P(t)~exp(-t/Ts). 
The calculations were performed for the emission times between and 200 fm/c and 
for freeze-out volume Vf = 3Vb, where Vq is the nuclear volume at normal density. 
The experimental filter was applied to be in line with the experimental definition of 
the correlation function. It is seen from Fig. 3 that the model calculations reproduce 
well the deep drop of the yield for small relative velocities. 

For the quantitative comparison of the data with the model calculations the 
Pirson criterion is used. The plot of and confidence level CL versus the mean 
decay time of the system is shown in Fig. 4. The mean life time of the fragmenting 
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Fig. 4. Comparison result of the measured and model calculated correlation functions versus 
emission time. CL - dashed line, - solid line. 
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system is equal to or less than 40 fm/c at the confidence level ^67%. This result is 
in agreement with the emis sion tim e values obtained in the previous studies of the 
thermal multifragmentatior{SEIHl9l_ 
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